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■ Abstract 

^ I In this article, we develop a new approach to functional quantization, which consists in discretizing 

O ■ only a finite subset of the Karhunen-Loeve coordinates of a continuous Gaussian semimartingale X. 

^0 ■ Using filtration enlargement techniques, we prove that the conditional distribution of X knowing its 

' first Karhunen-Loeve coordinates is a Gaussian semimartingale with respect to a bigger filtration. This 

allows us to define the partial quantization of a solution of a stochastic differential equation with respect 
to X by simply plugging the partial functional quantization of X in the SDE. 

Then we provide an upper bound of the L^'-partial quantization error for the solution of SDEs involving 
. the Z/^+^-partial quantization error for X, for e > 0. The a.s. convergence is also investigated. 

' Incidentally, we show that the conditional distribution of a Gaussian semimartingale X, knowing that 

(-H , it stands in some given Voronoi cell of its functional quantization, is a (non-Gaussian) semimartingale. 

As a consequence, the functional stratification method developed in [B] amounted, in the case of solutions 
of SDEs, to using the Euler scheme of these SDEs in each Voronoi cell. 
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J> ■ Cameron-Martin space, Wiener integral. 

00' 

'nJ- ■ Introduction 

in; 

Let (fi. A, P) be a probability space, and E a reflexive separable Banach space. The norm on E is denoted by 
I • I . The quantization of a i?- valued random variable X consists in its approximation by a random variable Y 
taking finitely many values. The resulting error of this discretization is measured by the norm of |X — y |. 
If we settle on a fixed maximum cardinal for Y(Cl), the minimization of the quantization error amounts to 
the minimization problem: 

• ■ 

^ ■ min {|||X - Y -.n^ E measurable, card(y(r2)) < n] . (1) 



A solution to ([T]) is an optimal quantizer of X. The corresponding quantization error is denoted by 
8n.p{X) := min {|||^ — ^lllp' Y -.Vt ^ E measurable, card(y(i7)) < A^j. One usually drops the p sub- 
script in the quadratic case (p = 2). This problem, initially investigated as a signal discretization method 
|10| . has then been introduced in numerical probability to devise cubature methods ^24j or to solve multidi- 
mensional stochastic control problems [3]. Since the early 2000 's, the infinite-dimensional setting has been 
extensively investigated from both constructive numerical and theoretical viewpoints with a special attention 
paid to functional quantization, especially in the quadratic case ^18j but also in some other Banach spaces 
|29j . Stochastic processes are viewed as random variables taking values in functional spaces. 
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We now assume that X is a bi-measurable stochastic process on [0, T] verifying E [l^tP] dt < +00, so 
that this can be viewed as a random variable valued in the separable Hilbert space L^([0, T]). We assume that 
its covariance function is continuous. In the seminal article on Gaussian functional quantization [18| . it is 
shown that in the centered Gaussian case, linear subspaces U of i^([0, T]) spanned by L^-optimal quantizers 
correspond to principal components oi X. In other words, they are spanned by the first eigenvectors of 
the covariance operator of X . Thus, the quadratic optimal quantization of X involves its Karhunen-Loeve 
eigensystem (e^, A^)„>i. If F is a quadratic A^-optimal quantizer of X and d'^ (N) is the dimension of the 
subspace of L^([0,T]) spanned by Y{n), the quadratic quantization error £'^{X) verifies 

£ff{X)^ Af +£2^ [(g)AA(0,Af) ] form>d^(iV). (2) 

£%{X) < Af + £jj ( (g) (0, Af ) j for 1 < TO < d^(iV). (3) 

j>m+l \j = l / 

To perform optimal quantization, the decomposition is first truncated at a fixed order m and then the 
R^-valued Gaussian vector, constituted of the m first coordinates of the process on its Karhunen-Loeve 
decomposition, is quantized. To reach optimality, we have to determine the optimal rank of truncation 
d-^ (N) (the quantization dimension) and the optimal (i'^(A'^)-dimensional quantizer corresponding to the 

first coordinates A/'(0,A^). A sharply optimized database of quantizers of univariate and multivariate 

Gaussian distributions is available on the web site www . quantize .maths-fi . com [26j for download. Usual 
examples of such processes are the standard Brownian motion on [0,T], the Brownian bridge on [0,r], 
Ornstein-Uhlenbeck processes and the fractional Brownian motion. In Figure [1] we display the quadratic 
optimal N- quantizer of the fractional Brownian motion on [0, 1] with Hurst exponent H = 0.25 and N = 20. 
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Figure 1: Quadratic iV-optimal quantizer of the fractional Brownian motion on [0, 1] with Hurst parameter 
H = 0.25 and N = 20. The quantization dimension is 3. 

From a constructive viewpoint, the numerical computation of the optimal quantization or the optimal 
product quantization requires a numerical evaluation of the Karhunen-Loeve eigenfunctions and eigenvalues, 
at least the very first terms. (As seen in [181 [TOl [^2j . under rather general conditions on its eigenvalues, the 
quantization dimension of a Gaussian process increases asymptotically as the logarithm of the size of the 
quantizer. Hence it is most likely that it is small. For instance, the quantization dimension of Brownian 
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motion with N = 10000 is 9.) The Karhunen-Loeve decompositions of several usual Gaussian processes have 
a closed-form expression. This is the case for standard Brownian motion, Brownian bridge and Ornstein- 
Uhlenbeck processes. The case of Ornstein-Uhlenbeck processes is derived in [61, in the general setting of 
an arbitrary initial variance tJo. Another example of explicit Karhunen-Loeve expansion is derived in [8j by 
Deheuvels and Martynov. 

In the general case, no closed-form expression is available for the Karhunen-Loeve expansion. For example, 
the K-L expansion of the fractional Brownian motion is not known. Yet, one can use numerical schemes 
to solve the correspnding eigenvalue problem. In [5j, the so-called "Nystrom method" is used to compute 
the first terms of the K-L decomposition of the fractional Brownian motion and to perform its optimal 
quantization. 

In this article, we propose a new functional quantization scheme for a bi-measurable Gaussian process 
X, which consists in discretizing a finite subset of its Karhunen-Loeve coordinates, instead of performing a 
full quantization. This partial functional quantization approach is motivated by two observations. The first 
one is that the conditional distribution of X knowing that it falls into a given Voronoi cell of its optimal 
quantizer is the crux of the recently developed functional stratification scheme |6j. It comes to conditioning 
the process with respect to its first Karhunen-Loeve coordinates. This work provides a better justification 
of the functional stratification scheme of [6]. The second observation is that one of the main purposes of 
the (full) functional quantization of X is to perform a quantization of the solution of a SDE with respect 
to X, when a stochastic integration with respect to X can be defined (see [25l [20l [27] ) ■ As (full) functional 
quantizers of X will typically have bounded variations, one needs to add a correction term to the SDE. 
Eventually, this comes to plug the functional quantizer of X in the SDE written in the Stratonovich sense. 
In contrast, the partial quantization of X can be directly plugged into the SDE written in the Ito sense. We 
provide a.s. and convergence results for this method. 

The paper is organized as follows: Section [1] provides background on quantization-based cubature formu- 
las which are needed for the following. In Section [51 we develop a notion of generalized bridge of a continuous 
Gaussian semimartingale which extends the generalized Brownian bridge introduced by Alili in jT]. We prove 
that under an additional hypothesis (fTZj) . the generalized bridge of a continuous Gaussian semimartingale 
remains a Gaussian semimartingale with respect to a bigger filtration and we derive its canonical decompo- 
sition. (Let us mention the thorough study of the properties of Gaussian semimartingales available in [13].) 
A similar result is stated when conditioning by a Voronoi quantizer. We pay a particular attention to the 
special case of generalized bridges that we call Karhunen-Loeve generalized bridges and which amounts to 
the conditioning of AT by a finite subset of its K-L coordinates. Section [3] is devoted to the partial func- 
tional quantization of continuous Gaussian semimartingales and its application to the partial quantization 
of solutions of SDEs. We finally give and a.s. convergence results for partially quantized SDEs. 

1 Quantization-based cubature and related inequalities 

The idea of quantization-based cubature method is to approach the probability distribution of the random 
variable X by the distribution of a quantizer Y of X. As y is a discrete random variable, we can write 

N 

Fy = J2 P^^v^ ■ If F : E ~^Ris a. Borel functional. 



Hence, the weighted discrete distribution {yi,Pi)i<i<N of Y allows one to compute the sum (|4]). We review 
here some error bounds which can be derived when approaching E[F(A))] by ([4]). See |i25j for detailed proofs. 

1. If X E , Y a quantizer of X of size N and F is Lipschitz continuous, then 



N 




(4) 



i=l 



\E[F{X)]-E[F{Y)]\ < [F]up\\X~Yh, 



(5) 
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where [-FJlip is the Lipschitz constant ofF. In particular, if (Yn)n>i is a sequence of quantizers such 

N 

that lim \\X — Yn\\2 = 0, then the distribution ^ pfS^i^ of Yn weakly converges to the distribution 
fx of X as N ^oo. 

This first error bound is a straightforward consequence of \F{X) ~ F{Y)\ < [F]up\X — Y\. 

2. If Y is a stationary quantizer of X, i.e. Y = E[X|y], and F is differentiate with an a-Holder 
differential DF for a S (0, 1], i.e. \DF{u) - DF{v)\l(e) < [DF]a\u - v\'\ for all {u,v) G E'^ where 
I • \l(e) is the operator norm on L{E), then 

|E[F(X)] - E[F(y)]| < [DFU\X - (6) 

In the case where F has a Lipschitz continuous derivative [a =^ 1), we have. [DF]i ~ [DF]i^ip. For 
example, if F is twice differentiable and D^F is bounded, then [Z3_F]Lip = ||-D'^-F||oo- 

This particular inequality comes from the Taylor expansion of F around X and the stationarity of Y . 

3. If F is a semi- continuou^convex functional andY is a stationary quantizer of X , 

E[i^(F)] < E[i^(X)]. (7) 

This inequality is a straightforward consequence of the stationarity property and Jensen's inequality. 

E[i^(y)] = E[i^(E[x|r])] < E[E[i^(x)|r]] = E[i^(x)]. 

2 Functional quantization and generalized bridges 
2.1 Generalized bridges 

Let {Xt)t<z[Q.T] be a continuous centered Gaussian semimartingale starting from on (f?,^, P) and its 
natural filtration. Fernique's theorem ensures that /q E dt < +00 (see Janson [14J). 
We aim here to compute the conditioning with respect to a finite family Zt ■= {Z^)i^i of Gaussian random 
variables, which are measurable with respect to a{Xt,t E [0,r]). (/ C N is a finite subset of N*.) As Alili 
in [T] we settle on the case where {Z!j.)i^j are the terminal values of processes of the form Zi = fi{s)dXs, 
i G I, for some given finite set / = {fi)i^i of Lf^^([0,T]) functions. The generalized bridge for {Xt)t£io,T] 

corresponding to / with end-point z = (zi)ie/ is the process {xf'^] that has the distribution 

V / te[o,T] 

X^-^ ^ C{X\Z'j,^ z.„ iel). (8) 

For example, in the case where X is a standard Brownian motion with |/| — I, f = {/} and / = 1, this is 
the Brownian bridge on [0, T]. li X is an Ornstein-Uhlenbeck process this is an Ornstein-Uhlenbeck bridge. 

Let H be the Gaussian Hilbert space spanned by {Xs)selo.T] and H-^^, the closed subspace of H spanned by 
(Z^)ig/. We denote by its orthogonal complement in H. Any Gaussian random variable G oi H can be 

orthogonally decomposed into G = Proj-g-^ (6*)+ Proj-^- (G), where Proj-g^^ and Proj-g- are the orthogonal 
projections on ifg-^ and iji . (Proji = /d// — Proj-g^^). With these notations, E [G| (Zy)ig/] = Proj;g^_^ (G). 
Other definitions of generalized bridges exist in the literature, see e.g. [23j. 

^In the infinite-dimensional case, convexity does not imply continuity. In infinite-dimensional Banach spaces, a semi- 
continuity hypothesis is required for Jensen's inequality. See |32| for more details. 
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2.2 The case of the Karhunen-Loeve basis 

As X is a continuous Gaussian process, it has a continuous covariance function (see [131 VIII. 3]). We denote 
by {ef , Xf)i>i its Karhunen-Loeve eigensystem. Thus, if we define function as the antiderivative of —ef 
that vanishes at t = T, i.e. f^{t) = ef{s)ds, an integration by parts yields 

X,ef{s)ds= I f^{s)dX,. (9) 



"'0 

In other words, with the notations of Section [^?Tl we have = Xgcf {s)ds =: Yi, the zth Karhunen-Loeve 
coordinate of X. 

For some finite subset / C N*, we denote by X^^y and call K-L generalized bridge the generalized bridge 
associated with functions {f^)iei and with end-point y — {yi)iei. This process has the distribution C{X\Yi = 
yi,i G I)- 

In this case, the Karhunen-Loeve expansion gives the decomposition 



X= Vr.ef + V ^Afe.ef, (10) 




where (^i)igN*\/ £^re independent standard Gaussian random variables. This gives us the projections Proj-g-^ 
and Proj^ defined in Section 12.11 It follows from ([TU]) that a K-L generalized bridge is centered on 
E = i e /] and has the covariance function 

r^l^(,s,i)-cov(X„X,)-^Afef(.)ef(0. (H) 
We have/i^r-^l^(t,t)(it J2 Af. 

Moreover, thanks to Decomposition PTI| . if X'-"^ is a K-L generalized bridge associated with X with terminal 
values y = {'yi)i^i, it has the same probability distribution as the process 



Y.y^ef{t) + x,~Y.( r 



X,ef{s)ds ef{t) 



This process is then the sum of a semimartingale and a non-adapted finite-variation process. 

Let us stress the fact that the second term in the left-hand side of ([TU)) is the corresponding K-L generalized 

bridge with end-point 0, i.e. Proj^ = X^-^. 

In [6], an algorithm is proposed to exactly simulate marginals of a K-L generalized bridge with a linear 
additional cost to a prior simulation of (^tp, • • • , ^t„), for some subdivision = to < ti ^ ■ ■ 'tn ~ T oi 
[0, T]. This was used for variance reduction issues. Note that the algorithm is easily extended to the case of 
(non-K-L) generalized bridges. 



2.3 GeneraUzed bridges as semimartingales 

For a random variable L, we denote by P [-jL] the conditional probability knowing L. We keep the notations 
and assumptions of previous sections. {X is a continuous Gaussian semimartingale starting from 0.) We 
consider a finite set / C {1,2, • • •} and {fi)iei a set of bounded measurable functions. Let X^^^- be the 
generalized bridge associated with X with end-point z — {zi)i^i. For i £ I, Z\ — fi{s)dXg and Zt = 
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Jirina's theorem ensures the existence of a transition kernel 

^zH((x.)..,„,,) :S(R^)xC"([0,.],R)^M+, 
corresponding to the conditional distribution £ {Zt\ ((^t)te[o.s])) • 

We now make the additional assumption fH|) that, for every s G [0, T) and for every {xu)ue[o.s] G ([0, s], M), 
the probability measure ^ZT|((Xt)te(o ]) i'^' i^u)ue[o,s]) is absolutely continuous with respect to the Lebesgue 
measure. We denote by njj,^)^^^ 7^ its density. The covariance matrix of this Gaussian distribution on 
writes 



Qis,T) := E [{Zt - E [ZT\{X^)uelo,s]]) {Zt - E [ZT\{Xu)ue[o,s]]y 



jij«e[o,s] 



If X is a martingale, we have (5(s, T) — yyj^ fiiu)fj{u)d{X)ujj^. -j^^a' We recall that a continuous centered 

semimartingale X is Gaussian if and only if {X) is deterministic (see e.g. ^28]). Hence, this additional 
hypothesis is equivalent to assume that 



Q{s,T) is invertible for every s G [0,T). 



The following theorem follows from the same approach as the homologous result in the article by Alili [Ij 
for the Brownian case. It is extended to the case of a continuous centered Gaussian semimartingale starting 
from 0. 

Theorem 2.1. Under the i FhI) hypothesis, for any s G [0,r), and for F-^^-almost every J G M"^, P ['IZt = 'z] 
is equivalent to P on J-^ and its Radon- Nikodym density is given by 

dP[-\ZT^z] ^ n(x„)„„„.,,,r(^) 

dr ~ Uo,t{z) 



Proof: Consider F a real bounded -measurable random variable and 
function. 



a bounded Borel 



• On the one hand, preconditioning by Zt yields 

E [F4>(Zt)] = E [E [F\Zt] (t)(ZT)] = J 0(z)E [F\Zt = z] Uo^Tiz)dz. 



(12) 



F / 0(z)n(x,).,[o,,j,T(^)d^ 



• On the other hand, as F is measurable with respect to J-^ , preconditioning with respect to J-^ yields 
E [F0(Zt)] = E [FE [(j)(ZT)\Tf]] = E 

Now, thanks to Fubini's theorem 

E [F^(Zt)] - J [Fn(^,)^^j„^,,^(z)] dz. (13) 



Identifying Equations p2p and (|13l) . we see that for P-g-^ -almost surely z G and for every real bounded 
-measurable random variable F, 



E [F\Zt ^z]=E 



F- 



n 



(14) 



no,T(^) 

Equation pi)) characterizes the Radon-Nikodym derivative of P \-\Zt — with respect to P on . □ 
We now can use classical filtration enlargement techniques [121 [T51 [5U] . 
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Proposition 2.2 (Generalized bridges as semimartingales) . Let us define the filtration Q-^'^ by Q^'^ := 
a (^Zx, J'^) , the enlargement of the filtration J-^ corresponding to the above conditioning. We consider the 

stochastic process D- :^ = ^^^^"^^ ^ ^ [O'^)' 

Under the fWf) hypothesis, and the assumption that is continuous, X is a continuous Q-^'^ -semimartingale 
on [0,T). 

Proof: is a strictly positive martingale on [0,r) which is uniformly integrable on every interval [0,i] C 
[0,r). Hence, as we assumed that it is continuous, we can write as an exponential martingale — 
exp {LI- \ (L-) J with LJ = ll {Df)-' dDI (as = 1). 

Now, as X is a continuous (J^'^, P)-semimartingale, we write X ~ + its canonical decomposition (under 
the filtration T^). 

• Thanks to Girsanov theorem, ;= M — (^M,L^^ is a (J^^,P \-\Zt = z] ) -martingale. 

— A consequence is that it is a [Q-^^^ .f" [■\Zt — z] ) -martingale. 

- And thus is a (^t/-^'^, P)-martingale. 

For more preciseness on this, we refer to P] Theorem 3] where the proof is based on the notion of 
decoupling measure. 

• Moreover, conditionally to Zt, V is still a finite-variation process V, and is adapted to Q^'^ . □ 

Remark (Continuous modification). In Provosition \2. ^ if one only assumes that has a continuous modi- 
fication!)^ , then with each one of its continuous modifications is associated a continuous Q^'^ -semimartingale 
on [0,T), and all these semimartingales are modifications of each other. 

Proposition 2.3 (Continuity of D^). If J-^ is a standard Brownian filtration, then D'- has a continuous 
modification. 

Proof: Consider s £ [0,T). Under the (j?Z|) hypothesis, the density n(x„)„£[o ^j,t writes 

^(x^)^e[o.s] = (S'T det Q(s, T))-^ exp {(z-E[Zt\ (Xu)uelo,s]] ) Qis, T)-'' (z -E[Zt\ {XuUelo.s]])*) ■ (15) 

Let us define the stochastic process H by Hg := E [ZT|(^M)iig[o,s]] ■ The so-defined process _ff is a J^-^ 
local martingale. Thanks the Brownian representation theorem, H has a Brownian representation and has 
a continuous modification. The continuity of s i~> det Q{s, T) and s i— >■ Q(s, T)~^ follows from the definition 
of Q{s, T) and the continuity of H (up to a modification). Hence, has a continuous modification. □ 

Remark. • The measurability assumption with respect to a Brownian filtration is satisfied in the cases 
of Brownian bridge and Ornstein- Uhlenbeck processes. 

• This hypothesis is not necessary so long as the continuity of the martingale Hg = E [-^tI (^u)«e[o.s]] 
can he proved by any means. 

2.3.1 On the canonical decomposition 

With the same notations, and under the (|W|) hypothesis, we can tackle the canonical decomposition of X^-'^. 
We have 
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and 



In n 



(x„)„„o,.,,t(^)) = -^ln(27rdetQ(s,r)) 



- i (z - E [Zt| (X„)„e[o,,]] ) Q(s, T)-' {z-E [ZT\{Xu)uelo.s]])* 



Using that for a positive continuous semimartingale S , d\n S — ^ (^^ ■ S") , we obtain 



dn 



.t(-) 



n 



(^u)ue(o,s] 



= din n 



finite- variation 
process 



= -id((z-E[ZT|(X„)„e[o,.]])Q(s,r)-iJz-E[ZT|(XO„e[o..]])*) + (f.-v. p.) 
- (dE [Zt| (X„)„e[o,.]] ) T)-i (z-E [Zt| (Xj„g[o,,]] )* + (f.-v. p.) . 

Hence, 

This expression can be further simphfied in the two following cases: 

• In the case where X is a martingale, owing to the definition oi Zj, we haveVj G /, E ^Zi^ iXu)ue[o.s] 
lo fj{'^)dXu so that 

d (X, L~)^ = (7(s)g(s, r)-i (z-E [Zt| (X,)„e[o,s]] )*) d{X)s 

- E /^(•^) E (z, -E [zf (x„)„e[o..]]) 



(16) 



As a consequence, Af-J^ ^ /,(s) ^ (Q(s,T)-i)^. (z, - E [z^ (X„)„e[o,,]] ) is a (g^>/,P [-I^t = z])- 

martingale. We have recovered Alili's result on the generalized Brownian bridge |T]. 

• In the case where the Gaussian semimartingale X is a Markov process, for every j (z I there exists 
gj £ L2([0,T]) such that E [z^| (X„)„e[o.,]] = fj{u)dXu + gj{s)Xs. Indeed, 

i-T 



E [Z^ (X„)„g[o,,]] = / fj{u)dXu + 

^ 



fj{u)dX^ 



Hence, if one assumes that {gj)j^j are finite- variation functions (which is the case when X is an 
Ornstein-Uhlenbeck process or a Brownian bridge), we have d ^X, E [Zt| (X,j)„g[o.]] )g = (/(■s) +5('S)) d{X)s, 
and thus 

d{X,L-)^ = ((7(.s)+5(.s))Q(,s,T)-i(z-E[z^|(X„)„g[o,.]l)*)rf(^). 

E (/.(s) +5.(5)) E U-'^Iz't (.XuUelo.s]]) d{X),. 



Example (Standard Brownian bridge). In the case where X = W is a standard Brownian motion with 
\I\ = 1, f = {/} and f = 1, Zt — Wt and W^'^ is a standard Brownian bridge. We have Q[s,T)^^ — 



and 



d{X,L^)t = ^ (z-E [WT\{Wu)uem]) 



dt = 



z-Wt 
T-t 



dt. 
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Thus, 



z-Wt Wt 
dWt = -dt + -dt + dWt. 

T-t T-t 



(g^'l .V[-\Wt =z]) -martingi 



ale 



^^k\ =0. Conversely, B e 
= Ik] - 0. Thus, 



The martingale part happens to be a [Q-^'^ ,^[-\Wt — z]j-standard Brownian motion, thanks to Levy's char- 
acterization of the Brownian motion. Thus we have retrieved the classical SDE of the Brownian bridge. 

2.3.2 Generalized bridges and functional stratification 

With the same notations, we set — Projp(Z7-) = ^ 7,1(7. (Zt) a stationary quantizer of Zt (where 

i=l 

r = {71, • • • , 7Ar} and C = {C'l, • • • , Cjv} are respectively the associated knots and Voronoi partition). 

Proposition 2.4 (Stratification). Under the FN]) hypothesis, for any s e [0,T), for any k £ {1, • • • ,N}, 
P ^Z^ = 7fe j > and the conditional probability P Z^ = 7/jj is equivalent to P on J-^ . 

Proof: Obviously, if A e is such that V[A] = 0, we have P 
J^f satisfies P = 7/jj = 0, then pre-conditioning by Zt, we get E 

V[B\'Zt = z] dP^^iz) = 0. Hence P[B\Zt = z] =0 for P^^ -almost every z £ Ck- 
Since P-g- (Cfe) > 0, there exists at least one element 1 £ Ck such that P [B\Zt = 2] =0. Now thanks to 
TheoremMll ^B] =0. □ 

Proposition 2.5 (Stratification). Let us define the filtration Q-^'^ by G^'^ '■— c {^t^ , the enlargement 
of J--^ corresponding to the conditioning with respect to Z^ . For k g {1, • • • , N}, we consider the stochastic 

process D]" ' j -^-^x for s £ [0,T). 

Under the IH]] hypothesis, and the assumption that D'"' is continuous, the conditional distribution 
C [x Z^^ of X knowing in which Voronoi cell Zt falls, is the probability distribution of a Q^'^ -semimartingale 
on [0,T). 

Proof: Using that P Z'" = "fkj is equivalent to P on , thanks to Proposition 12.41 we can mutatis 

mutandis use the same arguments as for Proposition 12.21 P [■\Zt — being replaced by P Z^ = 7fcj . 

D'^^ is a strictly positive martingale on [0,T) uniformly integrable on every [0,t] C [0,T). Hence, as 
0'**= is continuous by hypothesis, it is an exponential martingale D1'' = exp (ij*" — \ {LP"') ^, with Lj'' — 
lo i-^s'')^^ dD]'' (as Dq* = 1). Now, as X is a continuous (J^"^, P)-semimartingale, we write X = V + M 
its canonical decomposition (under the filtration J^'^). 

• Thanks to Girsanov theorem, AH" M - {M,L''"') is a (^T^ ,F [• Z^ = 7fe] ^martingale. As a 
consequence, it is a (^Q^'^,¥ ^- Z^ = 7fcj ^ -martingale and thus is a (tj^''", P) -martingale. 

• Moreover, conditionally to Z^ , V is still a finite- variation process V, and is adapted to Q^'^ . □ 

Proposition 2.6 (Continuity of D^''). If J-^ is a Brownian filtration, then D'^^ has a continuous modifi- 
cation. 

Proof: By definition, D'^'' is a J^^-local martingale on [0, T]. The conclusion is a straightforward consequence 
of the Brownian representation theorem. □ 

Considering the partition of i^([0,r]) corresponding to the Voronoi cells of a functional quantizer of X, 
the last two propositions show that the conditional distribution of the X in each Voronoi cell (strata) is 
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a Gaussian semimartingale with respect to its own filtration. This allows us to define the corresponding 
functional stratification of the solutions of stochastic differential equations driven by X. 

In [B] , an algorithm is proposed to simulate the conditional distribution of the marginals [Xt^ , • • • , Xt^ ) of 
X for a given subdivision = to < < • • • < in = T of [0, T] conditionally to a given Voronoi cell (strata) 
of a functional quantization of X. The simulation complexity has an additional linear complexity to an 
unconditioned simulation of {Xt„, ■ ■ ■ ,Xt^). We refer to [HI for more details. 

To deal with the solution of a SDE, it was proposed in [B] to simply plug these marginals in the Euler scheme 
of the SDE. Proposition 12.51 now shows that this amounts to simulate the Euler scheme of the SDE driven 
by the corresponding (non-Gaussian) semimartingale. 



2.4 About the (jH]) hypothesis 
2.4.1 The martingale case 

In the case where X is a continuous Gaussian martingale, the matrix Q{s,t) defined in Section 12.31 writes 

Q{s,t)={{rj.{u)Mu)d{x)^))^^^^^^^^. 

For 1 < s < t < T, the map {■]■) : {f,g) i— ^ J* f{u)g{u)d{X)u defines a scalar product on L^{[s,t], d{X)). 
Hence Q{s, t) is the Gram matrix of the vectors of i^([s, t], d{X)) defined by the restrictions to [s, t] of the 
functions {fi)ii£i. Thus, it is invertible if and only if these restrictions form a linearly independent family of 
L'^{[s,t],d{X)). (Another consequence, is that if Q{s,t) is invertible for some < s < t < T, then for every 
(w, u) such that [s,t\ C Q{u,v) is invertible). 

For instance, if X is a standard Brownian motion on [0,T], the functions {f^)i^j (associated with 
the Karhunen-Loeve decomposition) are trigonometric functions with strictly different frequencies. Hence, 
they form a linearly independent family of continuous functions on every nonempty interval [s,T) C [0,T). 
Moreover, the measure d{X) is proportional to the Lebesgue measure on [0, T] and thus Q{s, T) is invertible 
for any s G [0,r). Hence, the IH]] hypothesis is fulfilled in the case of K-L generalized bridges of the standard 
Brownian motion. 



2.4.2 Standard Brownian bridge and Ornstein-Uhlenbeck processes 



Brownian bridge and the Ornstein-Uhlenbeck process are not martingales. Hence, this criterion is not 
sufficient and the invertibility of matrix Q{s, T) has to be proved by other means. 
Following from the definitions of Q{s,T) and Zt, in the case of the K-L generalized bridge 



Q(s,T), 



E 



cov 



/J fi'{u)dXu - E [/J fi'{u)dXu {X^)uelo,s]] 

X (/J ffi^)dXu - E [/J /f (u)dX„|(X„)„e[o 
{Jjfi^{u)dxif\jjffiu)dxty), 



{Xu)ue[o,s] 



(17) 



where { Xu'' 



has the conditional distribution of X knowing {Xu)ue[o, 



s]- 



u<£[s,T] 

• When X is a standard Brownian bridge on [0,r], Xu"^ is a Brownian bridge on [s,T], starting from 
Xg and arriving at 0. 

It is the sum of an affine function and a standard centered Brownian bridge on [s, T]. 

• When X is a centered Ornstein-Uhlenbeck process, xi^'' is an Ornstein-Uhlenbeck process on [s, T] 
starting from Xs, with the same mean reversion parameter as X. 

It is also the sum of a deterministic function and an Ornstein-Uhlenbeck process starting from 0. 



10 



As a consequence, in these two cases, the quantity cov ^ {u)dXu \ fj^{u)dXu^^ can be computed by 

plugging either a centered Brownian bridge on [s, T] or an Ornstein-Uhlenbeck starting from instead of X^^^ 
in Equation (|17p. This means that Q{s,T) is the Gram matrix of the random variables (/^ f^{u)dGu) 

where the centered Gaussian process {Gu)ui^[s.T\ is either a standard Brownian bridge on [s, T] or an Ornstein- 
Uhlenbeck process starting from at s. Thus it is singular if and only if there exists {ai)i^i 7^ in such 
that 

l^ajf (u) ) dG„ = a.s.. (18) 



= -9M 

The case of Brownian bridge 

In the case where X is the standard Brownian bridge on [0,T], functions {f^)iei are C°° functions and 
G is a standard Brownian bridge on [s,T]. An integration by parts gives Gsg'{s)ds = a.s. and thus 
5' = on {s,T) and thus g is constant on [s,T]. The functions [f^)i^i form a linearly independent set of 
functions and, as they are trigonometric functions with different frequencies, they clearly don't span constant 
functions, so that Equation (jl8p yields ai — ■ ■ ■ — an = 0. Hence the IH]] hypothesis is fulfilled in the case 
of K-L generalized bridges of the standard Brownian bridge. 

The case of Ornstein-Uhlenbeck processes 

In the case where X is an Ornstein-Uhlenbeck process on [0, T], G is an Ornstein-Uhlenbeck process on [s, T] 
starting from 0. The injectivity property of the Wiener integral related to the Ornstein-Uhlenbeck process 

stated in Proposition 12.71 below, applied on [s,r], shows that Equation (|18p amounts to g ^ ([sJ"]-''*) q 
thus 

V- X L''(ls,T],dt) 

2_^a,f^ = 0. (19) 

Again, as {f^)i^i are linearly independent, we have ai = • • • = a„ = 0. Hence the IH]] hypothesis is fulfilled 
in the case of K-L generalized bridges of the Ornstein-Uhlenbeck processes. 

Proposition 2.7 (Injectivity of the Wiener integral related to centered Ornstein-Uhlenbeck processes). Let 
G be an Ornstein-Uhlenbeck process defined on [0,r] by the SDE 

dGt = -OGtdt + adWt with cr > and 6* > 0, 

where W is a standard Brownian motion and Go ^ Af (O, ctq) is independent of W . 
If g L^([0,T]), then we have 

fg{s)dG.^^ ^ g^'^t'^'o. 
Jo 

± /■* 

Proof: The solution of the Ornstein-Uhlenbeck SDE is Gt = Gqc^^* + ae^'-'^^UvVs. Hence, we have 

independent of VK 

g{s)dGs ^ -OGq [ g{s)e-^'ds+ j gis)dG°. 

Thus, by independence, if g{s)dGs = then g{s)dG'^ = 0. This means that we only have to prove the 
proposition in the case of an Ornstein-Uhlenbeck process starting from 0. 

We now assume that CTq = and we temporarily make the additional assumption that 6T < |. If 
g e L'^{[0,T]) and g{s)dGs = 0, then g{s)Gsds = a g(s)dWs, and thus, if T'^^ denotes the 
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covariance function of G, 

g{s)g{t)T"" {s,t)dsdt^(j'^ I g{sfds. (20) 



/Q "'0 

Applying Schwarz's inequality twice, we get 



[ g{s)g{t)T°^{s,t)dsdt< f g{sfdsJ f [ {VOU (^s,t)f dsdt. 
Jo Jo \ Jo Jo 



Hence, provided that 



^ ^ {ro^is,t))'dsdt<^, 



(21) 



Equality ^ implies g{s)^ds = 0. 

Now, we come to the proof of Inequality (|2ip . The covariance function of the Ornstein-Uhlenbeck process 
starting from writes 



For t € [0,r], we have (rO^^(s,t))' ds ^ ^ {2 ~ Ae-^'^'Ot - e^^eCT-t) _ 2e-2et _^ Se'^^^ - e-2''(^+*)) 
and thus 



-2BT ^ g~4eT\ 



^ ^ (rO^(.s,t))'dsd^=^(-5 + 40T + 80Te-2«^ + 4e 

Consequently, the function defined by 4){9) := (r'~'^{s,t)^ dsdt — writes 

1 

^ (-21 + 40r + SOer^'^^T + 4e"^^^ + e~^'^'^) . 



' 16 614 

Thus (j){e) < -16 + 126ir which leads to Inequality (gT]) thanks to the fact that OT < 

We now come back to the general case where we might have 9T > |. If this is the case, let us consider 
f := T - i, so that 9 (t - < |. For i e [f , t] , we have 

independent of {Ws) 



T 



e[T.T] =:G? 

The so-defined process (^t ) j,j ^ centered Ornstein-Uhenbeck process starting from and satisfying 
the same SDE as G. Hence, by independence, if /J" g{s)dGs = 0, then J~ g{s)dG'^ — 0. 

As (t - f ) < |, we can apply the result to (G?)^^ |-~ so that g^ ^ ^='^'' 0. If T6' < |, we then 

-L"^ ( [0 T] ) r '^1 

have g = 0. If it is not the case, we use the same method by using the decomposition of 0,r into 
0, T — I and \T — ^,T\ and so on. An easy induction finally shows that g 0. 

The converse is obvious. □ 

The case of a more general Gaussian semimartingale 

In Appendix |X1 we investigate the problem for more general Gaussian semimartingales. As we have seen in 
the case of the Ornstein-Uhlenbeck process, if functions {fi)^^i are linearly independent in L'^{[s^T],d{X)) 
for s S [0,T), the (fH) hypothesis comes to the injectivity of the Wiener integral with respect to X on 
span(/i)ig/ (on interval [s,T]). 
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3 K-L generalized bridges and partial functional quantization 



We keep the notations and assumptions of Section 12.21 As we have seen, Equation ([TU)) decomposes the 
process X as the sum of a hnear combination of the Karhunen-Loeve coordinates Y :— {Yi)i(zi and an 
independent remainder term. We now consider a stationary Voronoi iV-quantization of Y . Y^ can be 
written as a nearest neighbor projection of Y on a finite codebook T — (71, • • • , 7Ar). 

Y^ = Projp(y), where Projp is a nearest neighbor projection on T. 

For example, Y^ can be a stationary jjroduct quantization or an optimal quadratic quantization of Y . 
We now define the stochastic process X^-^ by replacing Y by Y^ in the decomposition ([TO)) . We denote 
X^-r ^ Proj,,p(X). 

iei ieN*\/ 

The conditional distribution of X^'^ given that Y falls in the Voronoi cell of 7^ is the probability distribution 
of the K-L generalized bridge with end-point 7^. In other words, we have quantized the Karhunen-Loeve 
coordinates of X corresponding to j S /, and not the other ones. 
The so-defined process X^'^ is called a partial functional quantization of X. 

3.1 Partial functional quantization of stochastic differential equations 

Let X be a continuous centered Gaussian semimartingale on [0,r] with Xq — 0. We consider the SDE 

dSt = b(t, St)dt + a(t, St)dXt, Sq = x M, and t e [0, T], (22) 

where h{t, x) and ij{t, x) are Borel functions, Lipschitz continuous with respect to x uniformly in t, a and 
6(-,0) are bounded. This SDE admits a unique strong solution S. 

The conditional distribution given that = for i g / of 5 is the strong solution of the stochastic 
differential equation dSt — b{t, St)dt + a{t, St)dXl'^, with Sq — x ^ R, and for t g [0, T], where X^"^ is the 
corresponding K-L generalized bridge. 

Under the (THI) hypothesis, this suggests to define the partial quantization of S from a partial quantization 
X^'^ of X by replacing X by X^'^ in the SDE (|22|). We define the partial quantization S^'^ as the process 
whose conditional distribution given that Y falls in the Voronoi cell of 7^,. is the strong solution of the same 
SDE where X is replaced by the K-L generalized bridge with end-point 7^. We write 

dSl'^ = b (t, Sl'^) dt + a [t, Sl'^) dX^'^. (23) 

Remark. The SDE is written in the ltd sens unlike in the previous works on full functional quantization 
J2T\ \ 25^ where the SDE was written in the Stratonovich sense. 

Here, the set / of quantized Karhunen-Loeve coordinates does not depend on the quantization level, while in 
the case of full functional quantization, optimality is reached by adapting the quantization dimension. The 
optimal quantization dimension (or critical dimension) has been thoroughly investigated in [T^, T9] and is 
shown to be asymptotically equivalent to the logarithm of the quantization level when in goes to infinity, in 
the cases of Brownian motion, Brownian bridge and Ornstein-Uhlenbeck processes. 

3.2 Convergence of partially quantized SDEs 

We start by stating some useful inequalities for the sequel. Then we recall the so-called Zador's theorem 
which will be used in the proof of the a.s. convergence of partially quantized SDEs. 
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Lemma 3.1 (Gronwall inequality for locally finite measures). Consider I an interval oj the form [a,h) or 
[a, 6] with a < b or [a,oo). Let fi be a locally finite measure on the Borel a-algebra of I. We consider u a 
measurable function defined on X such that for all t G X, J \u{s)\fi{ds) < +00. We assume that there exists 
a Borel function onT such that 



u{t) < Tpit) ■ 



/ u{s)^{ds), Vi e I. 



If 



either ijj is non-negative, 

or t ^ /i([a,t)) is continuous on I and for all t £l, \^jj{s)\fj,{ds) < 00, 
then u satisfies the Gronwall inequality. 

u(t) < ipit) + / exp(^([s, t)))n{ds). 

J[a,t) 

A proof of this result is available in fO", Appendix 5.1]. 

Lemma 3.2 (A Gronwall-like inequality in the non-decreasing case). Consider X an interval of the form 
[a, b) or [a, b] with a < b or [a, 00). Let fi be a locally finite measure on the Borel a-algebra of I. We consider 
u a measurable non-decreasing function defined on X such that for all t £ X, \u{s)\n{ds) < -\-oo. We 
assume that there exists a Borel function ijj on X, and two non-negative constants [A, B) S such that 



u{t) < i;{t) + A I u{s)n{ds)+BJ u{s)^ii{ds), \/teX. (24) 

[a,t) V ha-.t) 



If 



either tp is non-negative, 

or t ^ lJi{[a, t)) is continuous on X and for all t £ X, \i/j{s)\^{ds) < 00, 
then u satisfies the following Gronwall inequality. 



(t) < 2^P{t) + 2{2A + B^) I exp {{2A + B^) ^([s, t))) fi{ds). 

'a,t) 



Proof: Using that for {x, y) G M^, ,Jxy <\[^+ By), we have 



u{sf^Ji{ds) ] < I u{t) I u{s)fi{ds) ] <^ + ^ I u{s)fi{ds). 



a,t) / \ J[a,t) / 2i? Z 

Plugging this in Inequality (IMl) yields 

u{t) < 2ijj{t) + (2A + B'^) I u{s)^i{ds). 

Applying the regular Gronwall's inequality (Lemma 13. ip yields the announced result. □ 

Theorem 3.3 (Zador, Bucklew, Wise, Graf, Luschgy, Pages). Consider r > and X be aW'' -valued random 
variable such that X G L''"'"'' for some rj > 0. We denote by £N,r{X) the L^ optimal quantization error of 
level N for X, SnA^) := mm{\\X -Y\\^ ,\Y{n)\ < N}. 

1. (Sharp rate) Let Vx{dC) ~ 0(C)'^C + ^{dS,) be the Radon-Nikodym decomposition of the probability 
distribution of X . (v and the Lebesgue's measure are singular). Then if 

—+— 

£N,riX) - Jr,d>i{ (t)^{u)du) xN^i, 

where J^.d G (0,oo). 
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2. (Non-asymptotic upper bound) There exists Cd,r,i-i S (0, oo) such that, for every -valued random 
vector X , 

•iN > 1, SnAX) < Cd,rJ\X\\r+riN-i. 

The first statement of tlie tlieorem was first establisfied for probability distributions with compact support 
by Zador [31], and extended by Bucklew and Wise to general probability distributions on M'^ [4J. The first 
mathematically rigorous proof can be found in [11) . The proof of the second statement is available in [2T] . 
The real constant Jr^d corresponds to the case of the uniform probability distribution over the unit hypercube 
[0, ll-^. We have X,i = i(r + and J2,2 = (see tUj). 



3.2.1 LP convergence of partially quantized SDEs 



Lemma 3.4 (Generalized Minkowski inequality for locally finite measures). Consider X an interval of the 
form [a,b) or [a,b] with a < b or [a,oo). Let ^ be a locally finite measure on the Borel a-algebra of I. Then 
for any non-negative bi-measurable process X — (Xt)tei o.'^^d every p £ [l,oo), 

Xtti{dt) < [ \\Xt\\pti{dt). 
p Ji 

Proposition 3.5 (Burkholder-Davis-Gundy inequality). For every p g (0, oo), there exist two positive real 
constants Cp^^ and C^^'^ such that for every continuous local martingale {Xt)t£[o,T\ null at 0, 



„BDG 



< 



sup \Xs 

se[o,T] 



< C. 



BDG 



We refer to [25] for a detailed proof. 



Proposition 3.6 {L^ inequality). Let G be a standard Gaussian random variable valued in M. There exists 
a constant Gp > such that for every M > 1 



-MP-' exp (^-— j < E [\G\p1ig\>m] < CpM^-' cxp (^- — 



Consequently 



^ ) ^ Mi exp (-^) < ||GJ|G|>A/||p < (Cp)i/fMi exp (-^ 



where q is the conjugate exponent of p. 

Proposition 3.7 (The non-standard case and L^ reverse inequality). If H aG has a variance of , we 
obtain 



HI 



\H\>M\ 



< (J 



Gl 



|G|>^ 



^^(Cp)i/P(M)iexp 

^a-^Gpf/PMU^^i-^. 

\ Zva 



_M1 



(25) 



■■■VM 



Conversely, for some fixed rj > 0, and if M > 1, we have 



M > J-a^{p~l)W- 



qrj 



2q 



2 



(26) 



pa 



where yV_i is the secondary branch of the Lambert W function. For more details on the Lambert W function, 
we refer to f^. 
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Theorem 3.8 [L^ quantization of partially quantized SDEs). Let X he a continuous centered Gaussian 
martingale on [0,T] with = 0. Let S he the strong solution of the SDE 

dSt — b{t, St)dt + a{t, St)dXt, Sq = x, 

where h{t, x) and a{t, x) are Borel functions, Lipschitz continuous with respect to x uniformly in t, a and 
6(-,0) are hounded. 

We consider X^'^ a stationary partial functional quantization of X and S^'^ the corresponding partial func- 
tional quantization of S , i.e. the strong solutions of 

dSf'^ — b (j;, Sf'^^ dt + a (j;, S^'^^ dx('^ , Sq'^ = x. 

Then, for every p € (0,cx3), e > and t G [0,T), there exist three positive constants Cx,e,i, ^x.ej o,f^d 
Bx e I such that 



sup 

ve[o,t] 



Sy Sy 



<Cx.e,iexp\ ^X,e,/i - 





2q 




p+e 



B 



X,s,I 



Y -Y' 



p+s 



(27) 



where q is the conjugate exponent of p, where Y is defined from X hy Equation (iOj) and Y^ is the nearest 
neighbor projection on T. 

Remark. Using that yV_i(— a:) ^ lii(a:), we can see that the right-hand term in Equation 1127]) goes to 



as the quantization error 



Y -Y^ 



p+£ 



goes to 0. 



Proof: We decompose the process X into Xt ^ J2 Y,ef{t)+Xt' and X^-^ into Xt'' = J2 ^»^ef + 

where Y^ is the nearest neighbor projection of Y on F. For some fee {1, • • • , N}, conditionally to Y^ — jk, 
we have 

St-Sl'^^ f {b{u,Su)-b{u,Sl'^))du-^Y. f {'Jin, Su)- a {u,Si^^))YFdef{u) 

Jo jgj "'0 

+ f [Y,-Y[)a{u,Su)def{u)+ I {a{u, Su) - (t [u.S'J)) Gud{X)u 



This gives (conditionally to Y^ = jk) 



St-S, 



9 - '^^ ■ ^ 



«e[o,T] 



+ / [a{u,Su)-a[u,Si'^))dMy. 





du 



(ef)'(-) 


max 




)i: 


a _ al,r 











Hn,ax|/| max |(ef)'(w)|r^|y,-yr|+ f {a{u,Su)-a{u,Si'^))Gud{X)u 

riu,Su)-a[u,Si'^)))dMy 
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max 

i;6[0,t] 



IS 



< [6] Lip / max 

'q v£[0.u] 



As a consequence, conditionally to = jk, 

+ NLip|/| max (ef)'(w) 

ue[o,T] 

[aU.\I\ max | (ef )' h| |y, - 



max 
V iei 



max 

ve[o.u] 



du 



tte[o,T] 



max 

ve[o,t] 



cr(M, S'„) - a (u, S^''")) Gud{X)u 



max 



To shorten the notations, we denote, for a random variable V and a non-negligible event A, ||y||p,A '■= 
E[VP\A]^^''. Hence, using the Minkowski inequality and the generalized Minkowski inequality for locally 
finite measures (Lemma 13. 4p . we get 



max 



Sy Sy 



+ [a]Lip|/| max |(ef)'(u) 











/' 


max 


o _ oi,r 




lo 


velo,u] 





ue[o,T] 



[<j]up\I\ max (ef) {u)\t 
i^i I 



ie 

u€[0,T] 



iGl 



P,{Y^=lk} 



max 



du 

Sy Sjj' 



du 



max 

uG[0,t] 



max 




i>e[o,t] 





a(u,5„)-(T(w,5^'r))dM„ 



p,{yr 



Now, from the Burkholder-Davis-Gundy inequality, 



max 

ve[o,t] 



Sy Sy 



I,T 











f 


max 


Sy Sy ' 




'0 


i'G[0,u] 





du 



P,{Y^=lk} 



W]up\I\ max (ef)'(u) (maxr,^ ) 



u£lO,T] 



max 

velo.u] 



Sy Sy 



I,T 



du 



p,{y^=ik} 



[a]up\I\ max (ef) (u)|r 

iie[o,T] 



lei 



f 


a{u, Su) - (T ( 


Jo 





P:{y''=ik} 



p,{y=lk} 



c. 



BDG 



(o-(u, 5„) - CT (m, S'li''")) d{X)y 



(28) 



Now, from Schwarz's inequality 





< 


iei 





Y -Y' 



pAY^=lk} 



pdy^'yk} 



17 



(cr(u, Su) - O (^U,Si'^ j) Gy,l\G^\<M 

q _ ql.r 



From the generalized Minkowski inequality 

|G„|d(X)„ <j;; {a{u,S^)-a{u,Si^^))Gu 

= Jo - cr 5^'^)) G„l|G„|>j\/ + (ct(u, - cr G„1|g„|<m 

+ /o| 

< 2[CT]max/o 1 1 ^1 G„ I >M 1 1 ^ {yr^^^^ } d{X)u + M[a]Lip Jq 

We obtain, thanks to Proposition 13.71 

< 2[aU.^{X)t{Cpy/Pvf A'H exp (-^) +M[fT]Lip/o 
" V ^ 

where vf — max (Var(G„)). Moreover 

ue[o,t] 



d{X)y. 



d{X)u 
d{X)u 



d{X)u 



G _ qi,r 



pAy'^lk} 



d{X)u, 



<7{u,Su)-(j[u,Su^^^ d{X)u 



< 



pAY'-^Jk} 









f 


max 


q _ qlX 

Oy Oy 


Jo 


iei 


v£[0,u] 





d{X)u. 



P,{Y^='rk} 



Hence, Equation ([28)) becomes 



max 

vGlO,t] 



Sy Sy 



< MLip|/| max (ef ) (u) 

P,{Y^=-fk} ne[0,T] 



Y -Y' 



P,{Y^='(k} 



+ rjM 



[o-]Lip|/| max 
is/ 

ue[o,T] 



{ef)'{u) 



max 

ve[o.u] 

max 

\ lei 



q _ qi 

^y 



Y 



du 

Sy 5"^' 



max 

■ue[o,«" 



du 



G, 



BDG 



max 

iG/ 

i>e[o,ji] 



Sy Sy 



I.T 



P:{Y^=lk} 

1/2 



d{X)y 



pAY^=lk} 



+ M[<j]up max Sy - Si^^ d{X)y. (29) 
« 'Jo ^e[o."l P,{Y^=ik} 

We can then apply the "Gronwall-like" lemma 13.21 for locally finite measures to the non-decreasing function 

i/p 



sup 

v£[0,t\ 



Uy kjy 



pAY^=lk} 



E 


sup 


Sy Sy ' 


p 






_«e[o,t] 
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and with the locally finite measure fi defined by ii{du) = du + d{X)u, and we obtain 

i/p 



sup 

ue[o,t] 



Sv S^' 



< AfE 



< AfE 



Y -Y^ 



Y^=Jk 



Y' =7fe 
i/p 



+ mi) exp {{Ef'-'' + C^'*^) M([0,t))) 
+ 77m) exp (i?f'^''/i([0,t)))exp(C^'*V([0,t))) , 



where Ej is an affine function of niax|(7fc)j|. This yields 



sup 

ve[o,t] 



Sy 



i.r 



X 



E 



Y -Y' 



Y' 



Now, for e > and p = 1 + - and q = = 1 + 2 the conjugate exponent of p, we have, thanks to Holder's 
inequality 



E 



r r E 



Y -Y^ 



Y' 



< 



< 



E 



Y -Y' 



Y' 



E 





p+e' 







Hence, 



E 



< 



y^ 



E 



y _yi 



p+e' 



Now, as the so-defined function cf) is convex and as Y^ is a stationary quantizer of y, we have thanks to 

.p ■ 



Equation 0, \\<j){Y^) _ < ||0 (yf ||^ and 



'yl 



If one sets M = \ -Vt{p - l)W-i 



y_yr 



j'^yp j , where q is the conjugate exponent of p and >V_i is the 



secondary branch of the Lambert W function. Proposition 13.71 ensures that tim ^ t] 
finally have the following error bound 



sup 

ve[o,t] 



<Cx,e,iexp\ [a]Lip\l-Vt{p-l)W-l \ ; 

which is the desired inequality 







2q 




y _ yr 











„,,2^29/p 2g/p 
pVf Op Uj- 



y _ yi 



y _yi 



We 



p+e 



P+e 



□ 



Remark (Without the stationarity property). The last step of the demonstration of Theorem \3.8\ (the use 
of Jensen's inequality) relies on the stationarity of the quantizer Y. Now, without this stationarity hypothesis 
and under the additional assumption that 



rnB(o,i) 7^ 0, 



(A) 



we have for every i d I 



Y 



< 



\Y,\ < \Y\ 



Y^--f-° <2|y,|+7^ <2|y,| + l, where 7'''« e rn 5(0,1). 



Hence 



max 



Y, 



< 2max|yi| + 1. 

is/ 



19 



We notice that the function (t>{x) defined in the demonstration of Theorem 
some non- decreasing function ip. This implies 



writes 4'{x) = ip[ra.a,yiXi) for 



(j) (y) = tjj [m.a.yiYij < ij; (^max(2|yi| + 1) ) = 9!'(2|r| + 1). 



Hence, we can obtain the same conclusion as in Theorem[ 

Corollary 3.9 (i^ convergence). With the same notations and hypothesis as in Theorem \3.8\. consider 
/jsf^.r„j ^ sequence of partial functional quantizers of X and IS^'^"] the corresponding sequence of 
partial quantizers of S . (For n G N, r„ is assumed to have cardinal n.) 
If we make the additional assumption that the associated sequence of quantizers {y^" 
for the convergence for some e > 0, then for every t £ [0, T) we have 

i/p 

= O (n 



' nef 



is rate- optimal 



sup 

ue[o,t] 



Proof: As 



0, we have a.s. d [Y^" ^Y] 0. Hence, there exists Nq &N such that for 

every n > iVo, r„ satisfies Hypothesis From this observation, the result is straightforward consequence 
of the previous remark and Zador's theorem 13. 3[ which defines the optimal convergence rate of a sequence 
of quantizers. □ 



3.2.2 The a.s. convergence of partially quantized SDEs 

Theorem 3.10 (Almost sure convergence of partially quantized SDEs). Let X be a continuous centered 
Gaussian martingale on [0,T] with Xq — 0. Let S be the strong solution of the SDE 

dSt ^ b{t,St)dt + a{t,St)dXt, Sq ^ x, 

where b{t, x) and a{t, x) are Borel functions, Lipschitz continuous with respect to x uniformly in t, a and 
6(-,0) are bounded. 
We consider (X 



a sequence of partial functional quantizers of X and S ' " the corresponding partial 



functional quantization of S , i.e. the strong solutions of 



dS, 



b{t,St''^"]dt + a lt,S, 



dX. 



i.r„ 



(For n S N, r„ is assumed to have cardinal n.) We also assume that the sequence of partial quantizers of X 
is rate-optimal for some p > \I\, i.e. that there exists a constant C such that 



E 



Y -Y' 



i/p 



< Cn" 



for every 7i G N*, where Y is defined from X by Equation ilO\) and Y^ is the nearest neighbor projection on 
r. Then for every t £ [0,T), S^'^'^ converges almost surely to St. 

Proof: From Corollary [13 if t e [0,T), there exist r e (|/|,_p) and Nq eN such that for n > Nq, 











l/r 


E 


sup 




r 






«G[0,t] 









Hence, as |jy > 1, Beppo-Levi's theorem for series with non-negative terms implies 



E 



H sup 

„>i«e[o.t] 



Su S, 



< +00. 
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Thus 



E sup 

•i>i Me[o,t] 



< 



a.s. 



so that 



sup 

ue[o,t] 



-> P- 



a.s. 



□ 



Remark (Extension to semimartingales) . In Theorems \3.8\ and \3.1(A we limited ourselves to the case where 
X is a local martingale. The proofs are easily extended to the case of a semimartingale X as soon as there 
exists a locally finite measure v on [0,T] such that for every w G the finite-variation part dV{u!) in the 
canonical decomposition of X is absolutely continuous with respect to v. In particular, this is the case for 
the Brownian bridge and Ornstein-Uhlenbeck processes whose finite-variation parts are absolutely continuous 
with respect to the Lebesgue measure on [0, T] . 



A Injectivity properties of the Wiener integral 

In this appendix, we recall some resuUs on the definition of the Wiener integral with respect to a Gaussian 
process. We focus on the injectivity properties. Here, we pay special attention to the special case of the 
Ornstein-Uhlenbeck processes. 

The covariance operator and the Cameron-Martin space 

Consider X a bi-measurable centered Gaussian process on [0,r] such that ^[X'^]dt < oo and with a 

continuous covariance function T-^ on [0,r] x [0,T]. We denote by H span{Xf,te [0,T]} the 
Gaussian Hilbert space spanned by (^()tg[o,T]- The covariance operator Cx of X is defined by 

Cx: LH[0,T]) ^ L^{[0,T]) 

y ^ Cxy = ^[{y,X)X]. 

We have Cxy{t) =E[{y,X)X] (t) = E [j^ Xsy{s)dsXt\ = J^^ {t, s)y{s)ds where r^{t,s) = E[XtXs] is 
the covariance function of X. 

The Cameron-Martin space of X, (or reproducing Hilbert space of Cx), which we denote by Kx, is the 
subspace of i^([0, T]) defined by Kx ■= {t i~> E [ZXt] , Z E H}. Kx is equipped with the scalar product 
defined by 

{ki,k2)x ^E[ZiZ2] if fc, =E[Z,X], i = l,2, 

so that {Kx, {■)x) is a Hilbert space, isometric with the Hilbert space {{y,X) : y E L^{[0,T])} . Kx is 
spanned as a Hilbert space by {Cx(y) : y € L^([0,T])}. 

The Wiener integral 

Here, we follow the same steps as Lebovits and Levy-Vehel in [T7] and Jost in [TB] for the definition of a 
general Wiener integral. The difference here is that we use the quotient topology in order to define the 
Wiener integral in a more general setting. 

We define the map U : H Kx defined by U{Z){t) = E[ZXt]. By definition of H and Kx, C/ is a bijection 
and for any s E [0,T], we have U{Xs) = T'^{s,-). Consequently, Kx is spanned by (r'^(s, as a 

Hilbert space. Now, we linearly map the set of the piecewise constant functions £{[0,T]) to the Cameron- 
Martin space Kx by 

J: £{[0,T]) Kx 

i|Mi ^ r^{t,-)-T^{s,-), 

where |a, b\ stands either for the interval [a, b], (a, b), (a, 6] or [a, b). We equip £{[0, T]) with the bilinear form 
(•, •) J which is defined by 

{f,9)j--^ {Jf,J9)x- 
It is a bilinear symmetric positive- semirfe/iniie form. 

Remark. The so-called reproducing property shows that {l\o,t\, ^\o,s\ )j = r^(t, s) + r^(o, 0) - r^(o, s) - 
r^{0,t). When Xo = a.s., this gives j = {s,t). 
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Now, we define the equivalence relation ^ on £{[0, T]) hy x ^ y if {x — y, x — y) j = 0. On the quotient 

space E{[0, T]) := £{[0, T])/ ^, the bilinear form (•, •) j is positive-definite and thus it is a scalar product on 

E{[0,T]). In this context, J defines an (isometric) linear map from E{[0,T]) to Kx- Then, considering the 
completion F of E{[0,T]) associated with this scalar product, J is extended to F and U^^ o J : F ^ H is 
an (isometric) injective map that we call Wiener integral associated to X. 

T 

f{t)dXt -.^U'^oJif). 

Injectivity properties of the Wiener integral 

As we have just seen, the Wiener integral is an (isometric) injective map from F to H. Still, for example, 
when dealing with a standard Brownian bridge on [0,T], ||1[o^t]||j = 0, so that there are functions of 
f ([0,T]) which have a nonzero norm and a zero || • ||j norm. Injectivity only holds in the quotient space 
E{[0,T]) = £([0,T])/ - and its completion F. 

It is classical background that in the special case of a standard Brownian motion, || • ||,/ exactly coincides 
with the canonical norm so that F ~ i^([0,r]). 

Study of the case of Ornstein-Uhlenbeck processes 

From now, we will assume that X is a centered Ornstein-Uhlenbeck process defined on [0, T] by the SDE 

dXt = -eXtdt + adWt with ct > and 6* > 0, 

where is a standard Brownian motion and Xq ^ J\f (O, ctq) is independent of W. We make the additional 
assumption that 9T < |. The covariance function writes 

29 ' 

Proposition A.l (Semi-norm equivalence on £([0,r])). There exist two positive constants c and C such 
that for every f g £{[0,T]), c\\fh < \\f\\j < C\\fh. 

Proof: Let us consider / € £{[0,T]). We have 
11/11} =Var 1^-0^ f{s)Xsds + a f {s)dW^ < 2Y&T j ^ f{s)X,d^+2Y&i(^j^ f{s)dW, 

The solution of the Ornstein-Uhlenbeck SDE is Xt = Xge^^* + / cre^*''"*^^^/^- The so-defined process 

Jo 

V ' 

(Ar°)^gjQ is a centered Ornstein-Uhlenbeck process starting from 0. Hence we have 

11/11 2 < 2Var(Xo0/o^/(s)e-»^ds) -f 2Var(0/Q^/(s)XOds) -f 2Var(a/J'/(s)dW^,) 
< 202T Var(Xo) f{sfds + 2 Var [o fis)X^ds) + 2 Var {a fis)dW,) . 

As in the proof of Proposition \2.7l using that OT < 4/3, we can show that Var (6* f{s)X^ds) < 
Var(cr f{s)dWs). Hence 

\\m<{29^Tal+W) r fis fds, 
V 'Jq 
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which is the desired inequahty. Now we write 



i T / T \ T 

fis)dX,^-9 [ fis)Xoe-''ds+(-e [ fis)X°d.s)+a [ fis)dWs, 
"'O \ Jo J Jo 



= ■< =:G( 

where (Gq,G(,G2) is Gaussian and Gq is independent of G{ and Gj- Hence 



Var 



fis)dX,^ > Var [g( + gQ = Var [gQ + Var (g|) + 2cov (g{, G{] 



> Var (g{) + Var (g|) - 2y^Var (g{) Var (g{) = (^'Var (g{) - y'Var (g() j . (30) 

It has been shown at the beginning of the proof of Proposition 12.71 that there exists a constant K < 1 
independent of / such that Var (g() < KYar (g^). K was defined by 



JO 



where is the covariance function of the Ornstein-Uhlenbeck process starting from 0. Plugging this into 
Equation ((30)) yields 



Var 



^ fis)dX.^ > (l-\/i?)'var(G{) = (l-VK)\'\\f\\ 



This is the wanted inequality. □ 

A straightforward consequence of Proposition lA.ll is that ||/||j = ||/||2 — so that equivalence 
classes in £([0,T]) for the relation ~ are almost surely equal functions. Another consequence is that the 

sets of Cauchy sequences and convergent sequences for the two norms on E([0,T]) coincide, and thus the 
corresponding completions of E{[0,T]) are the same. In other words, in the case of Ornstein-Uhlenbeck 
processes that satisfy the condition 6T < |, we have F = L^([0,r]). 

The author is grateful to Benjamin Jourdain and Gilles Pages for their helpful remarks and comments. 
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